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Abstract: Inspired by the existing calculation oi B ^ ttK decays in supersymme- 
- try (SUSY), we evaluate the dominant SUSY contributions to 5° K^K~ . We 

■ show that the observables of this process can be significantly modified in the presence 

^ of SUSY. In particular, the branching ratio can be increased considerably compared 

to the prediction of the standard model (SM). The effect is even more dramatic for 
the CP- violating asymmetries Adir and Amix- These asymmetries, expected to be 
small in the SM {A^ir is predicted to take only positive values), change drastically 
with SUSY contributions. The measurement of these observables can therefore be 
used to detect the presence of physics beyond the SM, and put constraints on its 
parameters. 
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1. Introduction 

Current measurements of B decays show hints of physics beyond the standard model 
(SM), in CP- violating asymmetries in penguin-dominated b sqq transitions (g = 
u,d,s) [2], in triple-product asymmetries in i? ^ 07^*^, in the polarization 
measurements of B ^ 0-^*0 |, H and B pK* 0, ||, and in B ^ nK decays 



(branching ratios and CP asymmetries) |T0|, |TT], |12|, [13], |1J, |^ . These discrepancies 



are (almost) all not yet statistically significant, being in the l-2a range. However, 
if these hints are taken together, the statistical significance increases. Furthermore, 
they are intriguing since they all point to new physics (NP) in 6 — > s transitions. For 
this reason it is interesting to consider the effect of NP on B decays dominated by 
the quark-level h ^ s process. 

One such decay is 5° K^K~ . In the SM, its amplitude is given approximately 

by 

A{Bl ^ K+K-) = -P' -T' . (1.1) 

Here the prime on the amplitude stands for a strangeness-changing decay. In the 
above, P' and T' are the gluonic penguin amplitude and the color-favored tree am- 
plitude, respectively. These are estimated to obey the hierarchy P' : T' ~ 1 : A, 
where A ~ 0.2 [jl6|. There are other diagrams, but they are expected to be O(A^), 
and have been neglected above. 

The amplitude P' is actually composed of three pieces, P^, P^' and P/, where the 
subscript refers to the internal quark in the loop: 

p' = y:,yus k + k + y:,yts k 
^ v:^v,, (p; - p[) . (1.2) 
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In writing the second line, we have used the unitarity of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix to ehminate the V^lVts term, and we have dropped the 
V:,Vus term since jl/^K.! < IK^K.!- 

Note: even though the term V*,jVus {Pi — Pt) is at the level of other terms we 
have neglected, it can be retained by redefining the T' amplitude: 

T' ^ T' + KfeKs - P/) . (1.3) 

In the rest of the paper we will adopt this redefinition. Thus, T' has both a tree and 
a (small) penguin component. 

At the quark level, 5° K^K~ is described by 6 ^ suu. There are many poten- 
tial NP contributions, which at the quark level take the form [K^ K~\bT iS uT ju\B^g) , 
where the represent Lorentz structures, and color indices are suppressed. (We 
expect the size of all NP contributions to be at most of the order of \P'\-) This 
picture can be simplified by considering the strong phases. 

In Ref. it was observed that the NP strong phases are negligible compared 



to that of the (dominant) SM contribution P' . (Note: each NP contribution can in 
principle have a different strong phase.) Briefly, the argument goes as follows. All 
strong phases are due to rescattering from intermediate states, with a suppression 
factor of about 10-20. In the SM, the P'^ strong phase arises principally from the 
rescattering of the h ccs tree diagram, T^'. Since T^. is about 10-20 times bigger 
than P^, a strong phase of 0(1) is generated. By contrast, the NP strong phases can 
arise only from "self-rescattering," i.e. rescattering from NP operators themselves. 
As a consequence, these phases are only about 5-10% as large as that of P', and 
are therefore negligible. This leads to a great simplification: if one neglects the NP 
strong phases, one can combine all NP matrix elements into a single NP amplitude, 
with a single weak phase: 



E 



{K+K-\bTiSuTju\B'^^) = ^"e**" , (1.4) 



where the Qj are the coefficients of the operators and is the effective NP weak 
phase. 

Note that while this argument — that the NP strong phases are negligible — is 
quite general, there are still ways of evading this result. This can occur, for example, 
if certain NP amplitudes are larger than \P'\ and do not contribute to K~^K~ , 

but still contribute to the rescattering. This situation is perhaps unlikely, but the 
reader should be aware of these caveats. 

Note also that the T' strong phase is expected to be small. Thus, the relative 
strong phase between T' and the NP is small compared to that of P'. Below, we will 
take this to be (0 ± 10)°. 

In a previous article, three of us (DL, JM, JV) showed that one can measure 
the parameters |^"| and by combining measurements of 5° K^K^ and B^ 
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TT+TT" In the present paper, we consider the generation of \A^\ and within a 
specific NP modeh minimal supersymmetry (SUSY). 

Naively, one would guess that all NP contributions to \A^\ and are suppressed 
by M^/M'^p, where Mjvp ~ 1 TeV, and are therefore small. However, there are SUSY 
contributions involving squark-gluino loops. Since these involve the strong coupling 
constant a^, they are proportional to Og/M'^p, and so can compete with the SM 
contributions which are of order a/M^ [{as/a){M^/M'^p) ~ 1]. Thus, there are 
large SUSY contributions to the NP parameters. Indeed, these are the dominant 
effects, and are the only ones which are considered below. As we will see, one can 
generate an |^"| of the same order as \P'\, so that the amplitude for — * K^K~ 
can be written 

A{Bl -> K+K-) = -P' -T' + ^"e**" . (1.5) 

The effect of SUSY on the — > K^K'^ observables can therefore be sizeable, 
and we examine it here. We begin in Sec. 2 by establishing the SM predictions for the 
various observables in 5° K^K~ . In Sec. 3, we evaluate the SUSY contributions 
to the NP parameters \A^\ and (^u- With this information, in Sec. 4 we calculate 
the combined effect of the SM and SUSY on the 5° K^K~ observables. We note 
that the presence of SUSY can dramatically change the values of these observables. 
Thus, their measurements can both establish the presence of NP and constrain the 
SUSY parameter space. We conclude in Sec. 5. 



^0 ^ K+K-: SM Results 



s 



We begin with general definitions of CP-violating asymmetries. For the decay B^ — > 
/, where / is a CP eigenstate, one can measure two such asymmetries. The direct 
CP asymmetry takes the form 

Adir = o ! -,-2 1 (2-1) 
1^1 +1^1 

where A is the amplitude for B^ —>■ f. A is formed from A by changing the sign of 
the weak phases. The mixing-induced (indirect) CP asymmetry takes the form 

Im(e-^^^^*^) 

Amix — ~2 — -— 2 ' (2-2) 

1^1+ 1-^1 

where 0s is the phase of B^-B'^ mixing. 

We now turn to specific expectations for B^ — » K^K^ within the SM. This 
process has three observables: the two CP asymmetries mentioned above, and the 
branching ratio. Without calculation, we can estimate the expected size of the CP 
asymmetries. For 5° K^K~, since the amplitude T' is subdominant, to leading 
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order this decay is described by a single amplitude, V*^Vcs [P^ — -P/)- As such, in the 
SM the direct CP asymmetry is expected to be small, of order \T'/P'\ ~ A ~ 20%. 
Similarly, the mixing- induced CP asymmetry approximately measures 0^. Since 0^ is 
also expected to be very small (in the Wolfenstein parametrization ||19[, Im^<j ~ 5%), 
this asymmetry is expected to be correspondingly small. 

In order to calculate the SM predictions for these three observables, we need the 
magnitudes and relative weak and strong phases of P' and T'. The relative weak 
phase is 7, one of the three interior CP-violating angles of the unitarity triangle. 
This phase can be obtained from a fit to a variety of other measurements, some non- 
CP-violating. The latest analysis gives 7 = 

6ltr [H- Note that this error includes 



uncertainties in theoretical quantities. This value will be used in our analysis. 

For the magnitudes and relative strong phase of P' and T', we can proceed in one 
of two ways. One approach is to use a particular theoretical framework to calculate 
these quantities (see, for instance, [^, 0). Alternatively, one can use measurements 
of — ^ 7r+7r~, along with flavor SU(3) symmetry, to obtain P' and T' [O, |T^, E2|, E3 



vr+vr can be 



:2.3) 



In this paper, we adopt the latter approach. 

Neglecting small contributions, the amplitude for the decay 
written 

A{B',^n+n-) = -P-T . 
As above, we can write 

p = v:,Vud Pu + v:,v,a Pc + v:,Vtd Pt 

= K^Kd {Pu - Pt) + VXV,, (P, - Pt) . (2.4) 

The difference compared to P' is that one cannot neglect the first term. On the other 
hand, we can absorb it into the definition of T: 



T^T + K^Kd {Pu - Pt. 



(2.5) 



Thus, T is not a pure tree amplitude, but contains a penguin amplitude. 

As with S° K+K-, there are three measurements involving B^ tt+tt : 
the two CP asymmetries and the branching ratio. These suffice to determine the 
magnitudes and relative strong phase of P and T, given the knowledge of 7. Using 
flavor SU(3) symmetry, these can be related to the magnitudes and relative strong 



phase of P' and T' |12, 14, 22, 23 





r 




Vus 


n 




Y 






p'/r 




VcsVud 


P/T 




ycdVus 



:2.6) 



In the SU(3) limit we have TZc = 1, ^ = 1 and 9' = 6, where 6' and 6 are the relative 
strong phases of P' and T', and P and T, respectively. 
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Due to U-spin breaking, TZc gets both factorizable and non-factorizable contri- 



butions. The former have recently been calculated using QCD sum rules ||2J] and 
found to be sizeable: 

TZc = 1.76l°i? • (2.7) 

It should be noticed, however, that factorizable corrections are absent in the double 
ratio {P' /T')/{P/T). In our analysis, we use the central value of TZc- For the other 
quantities, we take ^ = 1.0 ± 0.2 (which we vary), and 9' — 9 = 0°. Whenever we 
refer to the U-spin limit, we will mean ^ = 1 and 9' = 9, but always taking the value 
of Eq. (|;7D for 7^c. 



With the experimental measurements of tt+tt^ and the theoretical values 

for the SU(3)-breaking parameters, we can obtain P' and T', which allow us to 
compute the SM expectations for the K^K~ observables. The latest B^ 

TT+TT" data is: 

(5.5 ±0.5) X 10"^ BaBar 
5i?(5° TT+TT-) = <j (4.4 ±0.7) X 10-^ Belle 

(5.0 ± 0.4) X 10"^ Average 

-0.09 ±0.16 BaBar [|7 
Adir{B^d ^^TT-) = { -0.52 ± 0.14 Belle [|| 

-0.33 ±0.11 Average 



( 0.30 ±0.17 BaBar ^| 
^mix(52 ^ TT+TT-) = \ 0.67 ±0.17 Belle ^| 

i 0.49 ±0.12 Average 

Regarding BRfJ^^, it is sometimes more useful to present the ratio of branching ratios 
of 5° ^ K+K^ and 5° ^ n+n--. i?^ = BR{Bl K+K')/BR{B^^ tt+tt-; 



The SM i?^ K^K~ predictions for all four quantities are shown in Table 1 (see also 
[0)- Obviously these values are correlated. Fig. |l] illustrates the main correlations 
between the observables, for different values of the SU(3) breaking parameter ^. 

Note that the CP asymmetries are allowed to take large values. This does not 
imply that our above argument about the expected smallness of these asymmetries 
is incorrect. Rather, it points to the largeness of the present experimental errors. 

Despite the large regions, there is still room for NP. If any of the correlations is 
found to lie outside of the allowed regions, this is a signal of physics beyond the SM. 

3. SUSY Contributions to and 

In this section we evaluate the SUSY contributions to \A^\ and We adopt the 
following procedure: 
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(xlO^) 


ps SM 


ASM 
dir KK 


ASM 
mix KK 


7 = 61° 

e = i 


(6.4,42.6) 


(1.2,9.3) 


(0.15,0.45) 


(-0.32,-0.10) 


7 = 61° 
^ = 1 ±0.2 


(4.2,61.9) 


(0.8,13.5) 


(0.12,0.56) 


(-0.38,-0.09) 


7 = (611D° 
e = l 


(5.0,60.7) 


(0.9,13.2) 


(0.08,0.58) 


(-0.34,0.08) 



Table 1: SM predictions for the branching ratio and mixing induced and direct CP- 
asymmetries. The impact of the uncertainty in the U-spin breaking parameter ^ and 
CKM-angle 7 is shown. 




Figure 1: Correlations between the observables A^^y — BR'^]^ and yl^" — Af^-^ 
for 7 = 61° and ^ = 1, ^ = 0.9 and ^ = 1.1. 

1. We consider all operators generated at the heavy scale, taken to be M^^. We 
compute the SUSY contributions to the coefficients of these operators. 

2. Using the renormalization group, we run the operator coefficients down to rrib. 
Operator mixing is included here. 

3. We compute the matrix elements of the various operators at m^,. This allows 
us to calculate |^"| and 

We closely follow the approach of Grossman, Neubert and Kagan (GNK) One 
difference is that GNK are interested in isospin- violating effects ("trojan penguins"), 
while we consider both isospin-conserving and isospin- violating contributions to \A^\ 
and Another difference is that GNK calculate the NP contributions to B ^ ttK, 
while we concentrate on K^K^ . Here the quark- level calculation is the same, 

and so our computation can be considered as a check. 

We begin by listing all the operators which are generated by the new physics at 
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the heavy scale. The NP effective hamiltonian is p 



TTNP 



V2 



.i,q=u,d 



c?(/i) Of + cf (/i) on + CM Qsa + CM Qsg 



(3.1) 



where 



01 
01 



{baSa)v-A{,qpqi3)v+A 
{baSa)v-Aiqf3qi3)v-A 
{ba(la)v-A{(l/3Si3)v+A 



01 
01 
01 



{baSf3)v-Ayqpqa)v+A , 

ipaSj3)v-A{,qiiqa)v-A , 

{haq(i)v-A{q(iS 

ajV+A 1 



(3.2) 



Qsg = (^?./87r2)m,6a^,(l - i^)G^^^s 



In the above, a and (3 are color indices, and the subscript V ± A indicates that the 
Lorentz structure between quarks is 7'^(1 ± 75). Despite the fact that, at the quark 
level, —>■ K^K' is 6 ^ suu, rf-quark operators must be included above since they 
mix with the u-quark operators upon renormalization to m^. Note that the above list 
includes the chromomagnetic operator Q^g. The operators Oj and Q^g are obtained 
from Of and Qsg by chirality flipping. 

The above list of operators includes new-physics contributions to electroweak- 
penguin operators. As we will see, these effects can be significant. This shows that, 
although the SM electroweak-penguin contributions to -B^ K^K^ are negligible, 
the same does not hold for the NP. 

We now must compute the NP contributions to the Wilson coefficients cf (/i) and 
cf(/i) at the heavy scale /i = Mw As discussed above, the dominant contribution 
comes from QCD penguin and box diagrams with squark-gluino loops. 



For our analysis, we follow Ref. ||3^ and take the general minimal supersymmetric 
standard model at the electroweak scale without assuming any flavor models at high 
energies (e.g. at the scale of grand unification). Here, the SUSY flavor-changing 
neutral current problem is avoided by assuming that the down squark is decoupled 
from the strange and bottom squarks That is, we write 

cosOrt -sine,e-"^4? 



dr 



L 



sin^ie*''^s 



+ cos 6'i.6° . 



(3.3) 



In the above, the superscript '0' indicates gauge eigenstates, and 5^ is a new CP- 
violating phase. There are similar expressions for the right-handed squarks. The 
\6l^r\ are taken to be < 45°. Similarly, the up squark is assumed to be decoupled 
from the charm squark, and up-top squark mixing can also be ignored [0. With 
these approximations, the Wilson coefficients c\q{^) and C5g(/i) vanish; the others 
are given in Appendix A. 

Once we have calculated, for given values of the SUSY parameters, the Wilson 
coefficients at My/-, the next step is to compute the renormalization-group running 
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of these, including operator mixing, down to mi,. The details of the computation are 
given in Appendix B. 

The final step in the program is to compute the hadronic matrix elements of 
the operators in Eq. ( p.2|) for K^K^ . These are calculated using the naive 

factorization approach. 

We define 

Al^{K-\{hu)x\Bl){K^\{us)y\Q) , (3.4) 

where X and Y refer to Lorentz structures. The pseudoscalar nature of the mesons 
implies that 



_ aV+A _ aV-A _ aV-A — A 
^V+A — ^V-A — ^V+A — ^V-A — ^ 



aS+P _ aS+P _ _aS-P _ _aS^P ^ q 
^S+P ~ ^S-P ~ ^S+P ~ ^s~p — ^ 



(3.5) 



which define the hadronic quantities A and S. After Fierz rearranging and factor- 
ization, the matrix elements of the operators read: 

(O3") = -vA , (02) = -A , (3.6) 
= A , {O^) = vA , 

where rj = I/Nq = 1/3. The matrix elements of the operators are just (O") = 

— {O^), the minus sign coming from the change A —A and S —S. Finally, we 
define 

X ^ -2S/A- ^ - (3.7) 

The NP amplitude can now be written as 



{K+K-\H-i\Bl) = ^ 



-x(^c^ + c?)-^(c?-4)-(c^-c^) 



A , (3.^ 



where the coefficients are evaluated at m^. The hadronic quantities x ^^id A can 
be calculated in terms of the meson masses, form factors and decay constants. Using 
the following expressions for the factorized amplitudes, 

iV2fKmj^ 



(i^+|M7^75s|0) = iV2fKPi^ {K+\u'y,s\0) 



22 1 



TTlu + rus 



qt nib 

(3.9) 

with Q- = ~ Qk- — 1k+ = P^y t^^^ 

X=— ^ -~1.18, 3.10 

A = iV2{ml - niDf-KF^-"^ ~ i 1.37 GeV^ . (3.11) 
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In the above, the values of the masses, decay constants and form factors are taken 
from Refs. |Tl|, |3ll . 



4. K+K -. SM + SUSY 

We are now ready to calculate the values of the various —>■ K^K^ observables in 
the presence of SUSY. We begin with BR^k and Adir kk- The parameters P' and 
T' are taken from the SM analysis (Sec. 2). The NP parameters \A^\ and are the 



modulus and argument of the amplitude in Eq. ( p.8|) . Finally, we must address the 
question of the relative strong phase of T' and A^. If the factor (P^ — P[) were not 
present in T' [Eq. (|1.3| )], we would say that the strong phase of T' is the same as 
that of the NP, i.e. it is negligible and 5^' — Sj^p = 0. However, (P^ — P/) is present. 
And since P^ can have a non-negligible strong phase due to rescattering from the 
b —>■ suu tree diagram, the relative strong phase of T' and can be nonzero. We 
take 5t' — S^p = (0 ± 10)°. The quantities BR^^ and Adir kk can now be obtained. 

The mixing-induced CP asymmetry, A^ix kk, can also be affected by the pres- 
ence of SUSY. However, in order to compute the allowed range, we must take into 
account the fact that this NP will also affect the P^-P^ mixing angle (ps- The SM 
predicts (ps ~ 0, because the combination of CKM matrix elements {V^*Vtb)'^ is real 
to a very good approximation |19|. On the other hand, in the SUSY scenario we 



consider, sizeable 0s is possible. Barring the simultaneous existence of LL and RR 
mixing, we get the following expression for 0^: 



arg 
11 



X 



75 -^Sig, 



(4.1) 



Here Xt = ml/m^ and Aj = Vt^V^*. The loop functions F and G, are given in 
Appendix A, and 5*0 is 

^ , , X* — 12x^ + 15x^ — 4:X + 6x^ Inx , , 

*»'^> = 4(131)3 (4-2) 

For the case of PP-mixing, we can use the same formula with L ^ R. This allows 
us to compute Amix kk in the presence of SUSY. 

The complete expressions for BR^k, Adir kk and Amix kk depend on a number of 
unknown SUSY parameters. These are the gluino and squark masses, and the angles 

and 5^,^^. (The relative strong phase 6t' — S^^p has been discussed above.) Our 
aim here is to see how the space of allowed values for the P^ — > K^K^ observables is 
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increased with respect to that of the SM alone (Sec. H). For this purpose, we take the 
following ranges/values for the SUSY parameters. For the angles, we take — 7r/4 < 
0L,R < 7r/4 and — vr < (5^ ^. < vr. For the masses we take rrig = ^ = ^ = 
250 GeV, 250 GeV < m^^^ < 1000 GeV, and 500 GeV < m,-^^ < 1000 GeV.' We 
also take m,^ = mg^. 

Some further constraints are imposed on the set of input parameters. First, the 
same SUSY contributions to —>■ K^K^ will also affect B nK decays. In 
particular, there will be effects on the quantities i?* and Acpi'ii'^K'^) l3^, whose 
definitions and values are ||32| 



^* = 2[BR{B^^7r^K^)+BR{B-^7r^K-)] " ^'^^ ^ ^^'^^ 

, , , BR(B+ TT+K°) - BRiB- -n-R^) , , 

Acpi-K^K^) = \ ^ = -0.020 ± 0.034 4.4 

' BR{B+ -K+K^) + BR{B- n-R^) ^ ' 



In order to incorporate these two constraints we follow the approach in Ref. pO | 
where QCD-factorization is used, except for the strong phase related to Aq^ which 
we take as a free parameter. Second, there are bounds from BR{B ^ + 7) and 



Am, [33 



2.92 X 10"^ < BR{B ^ X, + 7) < 4.12 x 10"^ (4.5) 
Am,/Amf ^ > 0.9797 (4.6) 



The measured values ( |4.3| )-( |^^ ) will therefore put additional constraints on the SUSY 



parameter space, which are taken into account in our analysis. In particular, the 
bounds (^]^) for BR{B — + 7) have a strong effect on the angle 6l, which for 
(^L,R = 7r/4 is restricted to the ranges 0.86 < Sl < 1.35 and 4.93 < 6l < 5.43. 

The NP amplitude is found to be small (even zero) for small w-squark masses 
{rriu ~ 250 GeV). However, for u- and s-squark masses close to 1 TeV and large s — b 
mixing, can be as large as 3.3 x 10"*^ GeV (see figure This number should 
be compared to the magnitude of the SM penguin amplitude \P'\ ~ 3 x 10^*^ GeV. 
Since these values are similar, SUSY contributions can yield important deviations 
from SM predictions. 

If we now allow for the variation of the SUSY parameters we find that Amix 
can take any possible value. Correspondingly, Adir can take any positive value, and 
negative values down to —0.5. A deviation from the SM range shown in Table 1 
could be explained within SUSY. 

Turning to the branching ratio, it can also receive a sizeable correction from the 
gluino contribution compared to the SM prediction. It can be almost 90% larger 
than the SM prediction in the U-spin limit. Even if one includes a large uncertainty 
of ± 20% from the U-spin breaking parameter ^, the super symmetric prediction for 

^We take a wider range for BR{B Xs + ^) to allow for various theoretical uncertainties. 
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300 400 500 600 700 800 900 1000 

rUu (GeV) 

Figure 2: versus the common mass ruu^ = niuj^, for mj^ ^ = m^^ ^ = rrig = 250 GeV 
and several values of mg^^, in the case of maximal s — b mixing {9r^l = 7r/4) and 
Sl-Sr = tt. 



BR%T (xl06) 


r)S susy 


A susy 
dir KK 


A susy 
mix KK 


(3.6,79.1) 


(0.7,17.2) 


(-0.5,1.0) 


(-1,1) 



Table 2: Allowed ranges of the B° ^ K+K' observables, including both SM + SUSY 
contributions. 

the — > X+X" branching ratio can be up to 30% larger than that of the SM in the 
same regions of the SUSY parameter space. The same applies to the ratio i?^*"*^. 

Table 2 summarizes all of these results. We see that there is a wide range 
in the values of the observables which are not allowed by SM but that are easily 
accommodated by minimal SUSY. We stress that this result holds even in a situation 
of quite constrained parameter space, large hadronic uncertainties and a ±20% of 
SU(3) breaking in ^. 

5. Conclusions 

At present, there are several hints of new physics (NP) in processes governed by 
6 — s> s transitions. For this reason, it is useful to consider the effect of NP on 
h s processes. One such decay is 5^ K^K~ . There arc many possible NP 
contributions to — > K^K~ decays. However, to a good approximation, all of these 
have strong phases which are small compared to those of the standard model (SM), 
and can therefore be neglected. In this limit, one can combine all NP contributions 
into a single term, parametrized by its magnitude and weak phase 
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In this paper, we have calculated the main supersymmetric (SUSY) contribu- 
tions to and $„, assuming that B'^ tt+tt" is unaffected. There are many SUSY 
effects. However, the principal ones come from squark-gluino loops, which involve 
strong couplings, and because (as/a)(M|,/M^p) ~ 1, they are not suppressed com- 
pared to the SM (Mjvp ~ 1 TeV). These are expected to be the dominant effects, 
and so we have included only these contributions. We have used naive factorization 
to compute the matrix elements and used data from S — > tttt decays to estimate the 
SM contribution. 

In the presence of such SUSY contributions, the predictions of the SM for — > 
K^K~ decays can be significantly modified, particularly for u- and s-squark masses 
close to 1 TeV and large s — h mixing. For example, we have found that the branching 
ratio can be increased. Even if one takes into account the large uncertainty due to 
the breaking of flavor SU(3) symmetry, the prediction of the SM + SUSY for the 
S° — > K^K~ branching ratio can be up to 30% larger than that of the SM alone. 

The situation is even more dramatic for the CP-violating asymmetries A^j^ and 
A^ix- In the SM, these are predicted to be small, with A^nr taking positive values 
only. On the other hand, in the presence of SUSY contributions, the range of Amix 
gets enlarged from —1 to 1, and Adir covers all the positive range, and also admits 
negative values forbidden to the SM. 

We therefore conclude that the study of i?^ K^K^ decays is very useful with 
respect to new physics. The measurement of its observables can be used to detect 
the presence of NP. Furthermore, the precise values of these quantities can be used 
to constrain the NP parameter space. In particular, this holds true for the case of 
SUSY, which can significantly modify the SM predictions. 
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Appendix A: Wilson Coefficients 

Using the approximations discussed in the text, the non- vanishing Wilson coefficients 
(WC's) at the heavy scale (M^^.) are given by 

1 5 12" 



ai sin 261^6"^^ 
4V2G^m? 



12 



7 1 3 2' 



6 



5 1 12 



'bLg 
sin2^^e*'^^ 



SLgJ 



17 3 



3 ^ OLg' 



where the functions F, G, A and -B are 

xlnx 



(5.1) 
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y\ny 



{x-y){x-iy {y~x){y-iy {x - l){y - 1) 
x^ In X In y 1 



+ 
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{x-y){x-iy {y-x){y-iy (x-l)(?/-l) 

1 (l + 2x)lnx 

2(1 -x)^ 6(1 -x)2 
ll-7x + 2x^ Inx 



18(1 -x)3 3(1 -x) 



(5.2) 



and Xg.g = m|/m?, where rug^ {q = d, u) is the mass of the i*^ squark mass eigenstate. 
The expressions for the coefficients are obtained from those in Eq. (|5.1| ) via the 
exchange L ^ R. Note that there is a relative sign difference between our c\ and 



that given in Ref. [^]. Our computation of the WC's agrees with Ref . [p^ . 

Using the same approximations, the SUSY contribution to the WC of the chro- 
momagnetic operator is given by 



r^susY, _ 

^t-^^'^fia — o . — [Is 



37r 



3 AV2Gpml 



hL-g) 



{hi 



sl 



(5.3) 



where the loop function is 

-11 + 51 X - 21 x^ - 19 x^ + 6 X (-1 + 9 x) log(x) 



SUSY, 



X 



72(-l + x) 



(5.4) 



Appendix B: Renormalizat ion-group evolution of the Wilson 
coefficients 



The QCD evolution of the Wilson coefficients (WC's) is given by 

C{^J) = U,{^i,M^)C{M^) (5.5) 
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where U^{fi, Myy) is the evolution matrix. Following GNK we work at leading order 
(LO)^ neglecting electromagnetic corrections to the anomalous dimension matrix of 
the operators. The chromomagnetic operator Q^g is included. The leading logarith- 
mic approximation depends only on the leading-order anomalous dimension matrix 



(0) 



35]. We perform the running translating the basis to the 12 AB = 1 SM 



7 



operators. 



10 



i=3 



(5.6) 

where Aj = V^^^s- In the above, and Q2 are tree operators, Qs-q are QCD 
penguin operators, and Qj-w are electroweak penguin operators: 



= (baUp)v-A{upSa)v-A Q2 

Q3 = (baSa)v-AT,q{qf3qf3)v-A Qi 

Qb = ibaSa)v-Aj2qiWp)v+A Q& 

Qt = \(baSa)v-AYliq^qi.W[^)v+A Qs 
Q9=^ ' 



{bo,Ua)v-A{upS(j)v-A 
{baSp)v-AY.q{Wc^)v~A 
{baS(3)v-A Y.qiWa)v+A 
libaSp)v-A J2q ^q{.Wa)v+A 



(5.7) 



YSa)v-A^q'^>liWp)v-A Qw = i(^aS/3)y_A^ 



)V-A 



75)G 



where Cg is the electric charge of quark q. 

In this basis the evolution matrix is 12 x 12, with the coefficients Cj related to 
the NP cf s in Eq. (|]1|) through 



'^l — 


-Xt (C5 + Cj) 




^1 — 


-At (C5 - 


Id) 




•-2 — 


-Xt (Ce + Cs) 




'-2 — 


-At (Ce - 


2*^8) 




•-3 ~ 


-Xt (C3 + Cg) 


^d 

•^3 ~ 


-At (C3 - 


2C9) 


C4 — 


—Xt (C4 + Clo) 


r-d 

C4 — 


-At (C4 - 


2*^10 



(5. 



Note that the Cg'g are zero at the scale in our case. We take them to be zero 



also at the scale, since the electroweak combination C5g(mb) 



-5,6 



'nih) is at LO 



a function only of c^'q{M^)^ and the QCD combination (c5g(mf,) + 2c5g(mfe))/3 is 
mostly dominated by the same combination at M^, taking into account that all NP 
penguin WC's are of similar size. 
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